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Abstract 

> 

. A photon of momentum k can have only two polarization states, not three. Equivalently, 

one can say that the magnetic vector potential A must be divergence free in the Coulomb 
gauge. These facts are normally taken into account in QED by introducing two polarization 
vectors e\(k) with A 6 {1,2}, which are orthogonal to the wave-vector k. These vectors must 
be very discontinuous functions of k and, consequently, their Fourier transforms have bad decay 
properties. Since these vectors have no physical significance there must be a way to eliminate 
them and their bad decay properties from the theory We propose such a way here. 

1 Introduction 

03 
S. 

In quantum electrodynamics it is necessary to choose a gauge for the electromagnetic field, and we 
. ^ | shall use the Coulomb gauge here for the reason that it is the only gauge in which one gets the cor- 

rect electromagnetic field and electromagnetic energy by a minimization principle. This comes from 
the fact that the magnetostatic interaction energy of a current distribution can be found by mini- 
mizing (87r) _1 J B 2 — J j-A with respect to the vector potential whereas the (positive) electrostatic 
interaction energy of a charge distribution p is not given by the minimum of 87r -1 f |V^| 2 — f 4>p. 
(For a discussion of this see, e.g., pQ.) In the Coulomb gauge the electrostatic part of the interaction 
among particles is given directly by Coulomb's law eiej/\xi — Xj\. The curl-free part of the electric 
field is not a dynamical variable in this gauge. 

The dynamical field, whose dynamics is 'quantized', is the magnetic field. The (ultraviolet 
cutoff) magnetic vector potential is customarily defined by 



A (*) = ^J2j R3 £ ~7W\ 2A(k) ( ax{kyk ' x + a ^ k >' ik ' x ) dk . M 



A=l ' 
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where the function Xa * s a radial function in k space that vanishes outside the ball whose radius is 
the ultraviolet cutoff A. The creation and annihilation operators of photons of momentum k and 
polarization A, a\(k) and aX(k), which act on Fock-space (over L 2 (M 3 ) (g>C 2 ), satisfy the canonical 
commutation relations 

[a\(k), a* u (q)} = 5(k - q)5 x ,u , [a x (k), a u (q)} = 0, etc . (1.2) 

The magnetic field is B(x) = curl A(x). The vectors £\(k) are two orthonormal polarization 
vectors, which are perpendicular to k as well as to each other. 
The field energy, Hf, sometimes called dT(u>), is given by 



Hf 



he V f \k\ a* x {k)a x {k)dk . (1.3) 



A=l,2 



There is no cutoff in Hf. The energy of a photon is /ic|/c|. 

The polarization vectors are necessarily discontinuous. They must be discontinuous on every 
sphere centered at the origin, in fact, for the well known reason that "one cannot comb the hair on 
a sphere" . One possible choice, but by no means the only one, is 



si(k) 



(fc 2 ,-fci,o) 



e 2 (k) = jj|||Aei(fc). (1.4) 

While polarization is physically meaningful and measurable, the polarization vectors and the 
corresponding operators aT(k) have no direct physical meaning. The polarization vectors merely 
form an arbitrarily chosen and hence unobservable basis for vectors perpendicular to k. It should 
be possible to define the theory without the unphysical operators aT[k). 

Why does this matter? For an atom in the ground state $ one expects that the photon density 
decays towards zero with increasing distance from the nucleus. Since all the relevant quantities 
of the radiation field are expressed in terms of operators that act on k space, it is convenient to 
establish this decay by finding a bound on 

2 r 

/ ||V*a A (fc)*|| 2 dfc , (1.5) 
A=l ^ 

noting that smoothness in k space translates into decay in configuration space. To calculate the 
above expression one needs to compute £\(k) whose singularity causes needless complications, 
as in 0E|. 

A related complication caused by the polarization vectors occurs if one wants to calculate the 
coupling function h\{x) that measures the strength of the interaction of the electron with the photon 
field. This function is given by the Fourier transform of the /c-dependent quantity appearing in 
1)1. 1J) . namely (with superscripts i = 1,2,3 denoting components) 



h\{y) = — j ^±=± e \{k)e-^ydk (1.6) 
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The reason for the interest in h\(y) is that in some problems, such as the verification of the 
binding condition for atoms [3] or the existence of the thermodynamic limit, it is necessary to 
localize the electromagnetic field in x— space. Thus, if we formally define 

ax(x) = — V -jr [ a x (k)e ikx dk (1.7) 
(2tH) J r3 

to be the Fourier transform of the operator a\(k) then 

2 

A\x) = W(x ~ ■) + d\{h\{x - ■) . (1.8) 

A=l 

(Formally, d\{h\{x — •) = J a\{y)h\{x — y)dy.) We want h\ to have a rapid fall-off in order that 
d\(y) be localized with y as close to x as possible. 

Unfortunately, because of the discontinuity of s\(k) it is very difficult to decide what the fall-off 
of h\{y) is. The decay, being slow in directions that are perpendicular to the direction of the 
singularity of s\(k), will be nonuniform. With a smooth cutoff function xa(&) we can ) by a suitable 
choice of the polarization vectors, make h\{y) decay in such a way that J \x\ 2 ^\h\(y)\ 2 dy is finite 
for all 7 < 1. This weak form of fall-off is useful, but inconvenient to work with. Morover, one can 
envision situations where the nonuniformity of the decay will cause problems. 

In contrast, if we omit the e\{k) then the fall-off of the basic coupling function 

will be 1 2/ 1 s/2 ag |^| _> ^ jf we choose the cutoff X\{k) to be smooth. To see this we note 
that |/c| -1 / 2 is the Fourier transform of |y| _5//2 in the sense of distributions Theorem 5.9]. The 
Fourier transform of xa is real analytic and decays faster than any inverse power of |y|. Hence, 
the convolution of xa (the Fourier transform of xa)> with |y| _5//2 decays like |y| -5 ^ 2 - As an aside 
we note that a sharp ultraviolet cutoff, h(y) would decay only like \y\~ 2 , which turns out to be 
insufficient for a good localization of the photon states. 

In an attempt to eliminate the polarization vectors from the formalism it was suggested in 
[3] that it would be better to start with a formalism that contains only "divergence-free" vector 
fields as the dynamical variables. In particular, the Fock space would be built over the L 2 -space of 
divergence- free vector fields instead of L 2 (M 3 ) <g> C 2 . While this approach removes the arbitrariness 
in the choice of polarizations it causes problems when one attempts to localize photons. If one 
multiplies a divergence-free vector field v by a smooth cutoff function j the result is not a divergence- 
free vector field in general. One possible localization procedure is to write v = curl w and then use 
the field jv + Vj A w. Again, this is fairly tedious. 

Since then, we have found an easier way and that is the subject of this paper. 

2 The Third Photon Mode 

Our proposal is really very simple. Let us introduce three photon degrees of freedom for each 
k E R 3 namely, 

a(jfc) = {^(k), a 2 (k), a 3 (k)}, (2.1) 
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with the o* satisfying canonical commutation relations [a l (k), aP^q)] = and [a 1 (A;), a J *(q)\ = 
SijS(k — q). Another way to say this is that we use the Fock space over L 2 (M 3 ) <g> C 3 . The field 
energy involves all three modes: 

3 

H f = h c y" / \k\ a 3 *(k)a J {k)dk . (2.2) 

The vector potential is replaced by 

A(a?) = — f — Lx A (A;)-^ A (a(k)e lk - x + a* (k)e' ik x ) dk , (2.3) 



2tt 

and the analogue of (jl.8|) is the vector equation 

A(x) = i curl k(h(x — •)) — i curl a*(h(x — •)) . (2.4) 
The function h is related to (|1.9j) 

r *> = h I p^"* ■ (2 ' 5) 

Note that 

i 3 

{curl a(h(x - ■)) } =^2{a j (d i h( y x--))-d i (d j h( y x--))^£ ije . (2.6) 

Following the argument about the fall-off of h(y), we see that the fall-off of h(y) is bounded 
by 1 2/ 1 3/2 an( j ^ g derivative are bounded by \y\~ 5 ^ 2 - Thus, the desired fall-off |y|~ 5 / 2 is obtained, 
using l|2.6jl . for the localization of A(x). 

We have the following situation: The quantized field operators have been increased from two 
to three for each k— value. Nevertheless, the vector potential A (which mediates the interaction 
of matter and radiation) still has the property that its k— Fourier transform is perpendicular to 
k. The field energy has been supplemented by an additional mode (which we might think of as 
'dark energy' since it does not interact with matter and is, therefore, not detectable). For reasons 
outlined at the beginning of this paper these additional modes are not longitudinal photons. Those 
have been eliminated from the theory by choosing the Coulomb gauge. 

What we shall show next is that this theory, in which the Hamiltonian describing matter and 
its interaction with radiation (relativistic or non-relativistic) is formally the same except for the 
extra invisible mode, gives the same physics as the old theory with the two polarization vectors. 

Before continuing, let us note that there is an alternative to (|2..3|) that accomplishes the same 
thing. We can replace the vector \k\~ 1 k A&(k) in (|2.3j) . whose ^-component is |/c| _1 A; 1 a^(k) e^, by 
the vector whose ^-component is 

E E 4 W 4 (*) o? (A:) = |E £a(*0 e x (k) ■ a(fc) J = a e (k) -^k- a(fc) . (2.7) 

The vector whose ^-component is displayed in (J2.7|) is the projection of a onto the plane perpen- 
dicular to the vector k. The vector |A:| _1 A; A a(/c) is a vector that is perpendicular to both a and k 
and whose norm is the same as the vector in (|2.7jl . 

The theories with (|2.7I) and with A a(fc) are indistinguishable. We shall continue with 

(|2~5|) . 
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3 Equivalence of the Two Theories 



Let us start with the theory defined by the vector potential ()2.3|) and field energy Hf in (|2,2|) . That 
is to say, the total Hamiltonian of matter plus radiation has the form 



-f^total — -ff matter (A) + Hf 



(3.1) 



where iJ matter (A) describes the matter (as particles, or as a quantized field, relativistic or non- 
relativistic) . The important point is that it depends on the radiation field only through the quan- 
tized vector field A in (|2.3[) . 

Instead of the three /c-dependent operators a we introduce 



a\(k) = e x (k) ■ a(k), 



k 



ao(k) = — ■ a(k) 



and check that (with v = A or v = 0) 

[a u (k), ct M (g)] = [a„(k), a* (q)] = S Vtfi S(k - q) . 



(3.2) 



(3.3) 



We also observe that A can be written using only the two a\(k) as in (jl.lj) and that the field 
energy (|2.2|) is 



Hf = hc / \k\ \ S~] a* x (k)ax{k) + ag(fc)a (fc) } dk 
Jr3 [t^i ) 



(3.4) 



All of this is breathtakingly elementary to verify. The conclusion, however, is interesting. 
This 'three-component' model (|3.1[) . which is simpler to deal with than the usual 'two-component' 
model because the polarization vectors are absent, has the property of merely describing the 'two- 
component' theory plus one totally independent scalar field ao whose time evolution is governed by 
the Hamiltonian 

H = he I \k\a* (k)a {k) dk . (3.5) 

Thus, the (Heisenberg) time evolution of ao is simply ao{k,t) = ao(/c)e lc ' fc '*. 

The eigenvalues of (|3.1|) . on the other hand, are those of the original 'two-component' model 
plus the scalar field energy, whose eigenvalues (if the radiation field is enclosed in a finite box) are 
of the form 

he \k\ n k ■ (3-6) 

The are nonnegative integers, of which only a finite number are positive. For the ground state 
we choose all = 0. The energy in the scalar field mode is not observable. 

If we are interested in the thermodynamic limit of ordinary matter coupled to the radiation field 
in a positive temperature Gibbs state we have to proceed as follows. First, we imagine the universe 
to be a huge box of volume V, while the matter is confined to a much smaller box of volume V. We 
would have to imagine this, even if we stayed with the usual 'two-component' formulation. Then 
we would take the limit V — > oo, after subtracting the well known positive temperature free energy 
F(T) of the field, whose asymptotic value was calculated by Planck in 1900 jH]. 



F(T) 



^Bolt z 



■ v 

'(2vr) s 



log ( 1 — exp 



hc\k\ 



^Bolt z 



dk 



(3.7) 
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This subtraction is necessary in either theory, the only difference being the subtraction of the scalar 
field energy in our 'three-component' theory. Hence the factor 3 instead of 2 in (|3.7|) . After this 
V — > oo limit, one takes the usual V — ► oo limit. In the end, the scalar field contributes nothing. 
Its role is only to contribute some simplification to a difficult calculation. 
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